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ABSTRACT 

Combining resolved optical spectroscopy with panoramic HST imaging, we study the dynamical properties 
of spiral galaxies as a function of position across two intermediate redshift clusters, and we compare the cluster 
population to field galaxies in the same redshift range. By modelling the observed rotation curves, we derive 
maximal rotation velocities for 40 cluster spirals and 37 field spirals, yielding one of the largest matched 
samples of cluster and field spirals at intermediate redshift. We construct the Tully-Fisher relation in both V 
and K$ bands, and find that the cluster Tully-Fisher relation exhibits significantly higher scatter than the field 
relation, in both V and K$ bands. Under the assumption that this increased scatter is due to an interaction 
with the cluster environment, we examine several dynamical quantities (dynamical mass, mass-to-light ratio, 
and central mass density) as a function of cluster environment. We find that the central mass densities of star- 
forming spirals exhibit a sharp break near the cluster Virial radius, with spirals in the cluster outskirts exhibiting 
significantly lower densities. We argue that the lower-density spirals in the cluster outskirts, combined with the 
high scatter in both K$- and V-band TF relations, demonstrate that cluster spirals are kinematically disturbed 
by their environment, even as far as 2Rym from the cluster center. We propose that such disturbances may be 
due to a combination of galaxy merging and harassment. 

Subject headings: galaxies: clusters: individual (CI 0024+1654, MS 0451-03) — galaxies: spiral — galaxies: 
evolution — galaxies: kinematics and dynamics — galaxies: stellar content 



1. INTRODUCTION 

The observed tight correlation between the rotation veloc- 
itie s of spiral galaxies an d their total luminosities, first noted 
by iTullv & Fisherl (fl977h . has proven invaluable in helping 
to pin dow n the extragalactic di stance scale in the local uni- 
verse (e.g. ITullv & Piercel 120001) . Since then, many authors 
have attempted to leverage the so-called Tully Fisher relation 
to study the evolution of spiral galaxies as a function of red- 
shift, generally interpreting any deviation from the local Tully 
Fisher (TF) relation as an evolution in luminosity. 

These studies have yielded mixed results, however, with 
conflicting estimates of the rate of B-band evolution as a 
function of redshift dBamford et all 120061: iBohm et alll2004t 
iMilvang- Jensen et alJl2003b IVogt et al.ll 19961) . Recently, sev- 
eral authors have documented many of the systematic er- 
rors that make comparisons between studies very diffi- 
cult dMetevier et all 120061: iNakamura et all 120061) . Indeed, 
iNakamura et all (l2006h argue that the most certain method of 
using the TF as a measure of evolution is to construct a large 
matched sample of galaxies, consisting of nearly equal num- 
bers of cluster and field galaxies all measured in the same 
way. 

In a similar manner, large matched sa mples such as 
those p resented by Nakamura et al. (2006) and Bamfor det all 
(2005) can also be effectively used to measure a different sort 
of spiral galaxy evolution: that caused by infall into a galaxy 
cluster. By carefully selecting galaxies across a wide range of 
environments, in and out of clusters, we can gain a better un- 
derstanding of the changes in star formation and kinematics 
that a spiral may undergo as it falls into a cluster potential. 
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In this paper, we attempt to construct such a sample out of 
a large survey we are performing of two massive galaxy clus- 
ters at intermediate redshift: CI 0024+17 (hereafter CI 0024) 
at z= 0.39 and MS 0451-03 (hereafter MS 0451) at z = 0.54, 
which were selected to be complementary in their global 
properties, as part of a larger project to understand the role 
of the cluster environment in galaxy evolution. We make use 
of high quality spiral rotation curves determined from Keck 
spectroscopy to measure the Tully Fisher relation. The sam- 
ple is large enough to allow a first investigation of the scatter 
of the relation as a function of cluster-centric radius. We also 
study a control sample of field galaxies in a range of redshifts 
centered about the cluster redshifts, in order to asses differ- 
ences between field and cluster spirals. Our sample of 40 
cluster galaxies is the largest yet reported with an associated 
field sample (37 galaxies), and provides a powerful means to 
examine the effect of the cluster environment on infalling star- 
forming spirals. 

In order to disentangle environmental processes affecting 
the dynamics of infalling spirals and those affecting their stel- 
lar populations, we contrast the trends in the TF relation with 
those observed for integrated V—K$ colors and disk mass den- 
sity. The latter quantity is constructed from rotation curves 
and HST scale lengths, and is predic ted to be sensitive to the 
strength of the 'harassment' process dMooreet al.lll999l) . 

A plan of the paper follows. In § 2, we describe our clus- 
ter survey and sample selection, and then detail our procedure 
for deriving maximal rotation velocities in § 3. In § 4, we 
present our results on the cluster and field Tully Fisher rela- 
tion, with an examination of other dynamical quantities in § 5. 
The Tully Fisher rela tion in CI 0024 has also been studied by 
iMetevier et all (2006), and in § 4 we also directly compare ro- 
tation measurements for several galaxies in common between 
studies. In § 6 we discuss these results in light of proposed 
physical mechanisms acting in the cluster environment. In 
this paper, we adopt a cosmology with f2 m = 0.3, £!a = 0.7, 
and H Q = 70km s^Mpc -1 . 
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2. OBSERVATIONS 

In this study, we leverage a large imaging and spectroscopic 
survey of two massive intermediate redshift clusters: CI 0024 
at z = 0.39 and MS 0451 at z = 0.54, yielding high quality 
spectra of both field and cluster spirals, suitable for extracting 
rotation curves. In the following, we describe our photometric 
and spectroscopic observations, data reduction, and sample 
selection. 

2.1. Imaging 

We make use of HST imaging of CI 0024 and MS 0451 from 
the co mprehensive wide-field survey described in iTreu et alj 
(2003) and Smith et al. (2006, in preparation). In CI 0024, 
HST coverage consists of a sparsely-sampled mosaic of 39 
WFPC2 images taken in the F814W filter (~ I band), provid- 
ing coverage to a projected radius > 5 Mpc. MS 0451 ob- 
servations were taken with the ACS, also in F814W, and pro- 
vide contiguous coverage within a lOMpcx lOMpc box cen- 
tered on the cluster. Both sets of observations are complete to 
FSUW> 25. 

For CI 0024, ITreu et al.l (120031) reported reliable morpho- 
logical classifications to 7 = 21.1. The MS 0451 observa- 
tions are proportionately deeper, and so reliable morpholog- 
ical classification is possible to the same rest frame absolute 
magnitude (My = -19.5). All g alaxies brighter th an this limit 
are classified visually following lTreu et all d2003l) . 

We also use panoramic ground-based /G-band imaging of 
both clus ters, and J-band ima ging of CI 0024, with the WIRC 
camera dWilson et al.l 120031) on the Hale 200" Telescope. 
These data comprise a 3 x 3 mosaic of pointings, spanning a 
contiguous area of 26' x 26' centered on each cluster. Obser- 
vations were made in 2004 November for MS 045 1 and 2002 
October for CI 0024. The d etails of the obser vations and data 
reduction are described by Smith et al. (2005) and Smith et 
al. (2006, in prep.) for CI 0024 and MS 0451, respectively. 
Point sources in the final reduced mosaics have a full width 
half maximum of 0.9" and 1.0" in CI 0024 and MS 0451 re- 
spectively, and the 3 — a point source detection thresholds are 
J = 21 .6 and K s =19.7 for CI 0024 and K s = 20.2 for MS 045 1 . 

These near-infrared data are supplemented with wide-field 
ground-based optical imaging. We make use of BVRI-band 
imaging of CI 0024 with the 3.6-m Canada-France-Hawai i 
Telescope using the CFH12k camera (Cuillandre et al.ll2000l) . 
full details of wh ich are available in lCzoske et al.l d2002l) and 
ITreu et al l (120031) . MS 0451 was observed bv lKodamaet al l 
d2005l) for the PISCES survey through the BRI-band filters 
using Suprime-Cam on the Subaru 8-m Telescope. Fu ll de- 
tails of these data are published by Kod ama et al.l (120051) . The 
CFH12k data reach 3-cr depths of B = 27.8, V = 26.9, R = 26.6 
and 1 = 25.9 in ~ 0"9 seeing, and the Suprime-Cam data reach 
3-er depths of B = 28.1, R = 27.3, 1 = 25.8 in seeing ranging 
from 0."6 to 1". The field of view of all of these optical data 
is well matched to the area surveyed by the near-infrared mo- 
saics discussed above. 

2.2. Spectroscopy 

Observations with the DEIMOS spectrograph on Keck II 
from October 2001 to October 2005 secured spectra for 
over 500 members of both CI 0024 (0.373 < z < 0.402) 
an d MS 045 1 (0.52 < z < 0.560). D etails are provided 
in ITreu etal.1 (120031) . iMoran et al] (120051) . and Moran et al. 
(2007, in prep). Briefly, we observe with 1" x 8" slits, with 
a typical velocity resolution of 50 km s" 1 . For each cluster, 
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FIG. 1. — Top: Postage stamp image of a 2D spectrum, centered on the 
observed wavelength of Ha for the galaxy N57426, in the field of MS 045 1 . 
The solid line across the image indicates the spectral row that is plotted in 
the lower panel. Bottom: One spectral row from the 2D spectrum. The solid 
curve indicates the Gaussian fit to this row. The stellar continuum emission 
has been subtracted according to the method described in the text. 

spectral setups were chosen to span rest frame wavelengths 
from ~ 3500A to ~ 6700A, covering optical emission lines 
[OII],[OIII], H/3, and, more rarely, Ha. Exposure times to- 
taled 2.5hrs in CI 0024 and 4hrs in MS 0451. In this spec- 
troscopic survey, we have also obtained spectra for over 2500 
field objects, with 700 having redshifts similar to the clusters 
(0.3 <z< 0.65). 

DEIMOS data were r educed using the D EEP2 DEIMOS 
data reduction pipeline dDavis et alj 120031) . which produce 
sky-subtracted, wavelength-calibrated two-dimensional spec- 
tra, suitable for identifying and extracting extended optical 
emission lines. Redshifts for all galaxies were verified by eye. 

2.3. Sample Selection 

Our current sample is drawn from the set of surveyed galax- 
ies with both HST imaging and available spectra. From the 
HST imaging, we select candidate galaxies that are morpho- 
logically classified as spirals (T-types 3, 4, or 5). For ease 
of comparing cluster spirals to field, we construct a matched 
sample of field spirals, all with HST imaging from the CI 0024 
or MS 045 1 mosaics, and selected to lie in a redshift range that 
brackets the two clusters, 0.3 < z < 0.65. 

In our spectroscopic survey, targets were selected randomly 
from an F814W-limited sample, to FSUW = 22.5 in the field 
of CI 0024, and to F814W = 23.0 in MS 0451. Here, we 
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FIG. 2. — Postage stamp images and rotation curves for cluster spirals, arranged in order of ascending projected radius, R, normalized by each cluster's Virial 
radius Rym- HST images are 9"", in F814W. Observed rotation curves are plotted as error-barred points, solid lines represent the best-fit rotation curve model, 
and the quality code is displayed on the plot. Radii in kpc are measured along the semimajor axis. In case where more than one emission line yielded a rotation 
curve, we plot the best one. On the postage stamp images, the position of the 1" X 8" slit is indicated by the solid black box. Some slits are not aligned with the 
galaxy major axis. At the redshift of the clusters, the spectroscopic seeing of 0"8 is equivalent to ~ A.5kpc. 
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Fig. 2. — (Cont.) 
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Fig. 3. — Postage stamp images and rotation curves for all field spirals, arranged in ascending redshift order. 



focus on a subset of these galaxies that have been observed 
with slits aligned along the galaxy major axis, in order to se- 
cure resolved spectra of spirals with extended emission lines. 
The CI 0024 campaign primarily targeted known cluster spi- 
rals for spectroscopy along the major axis, but serendipitous 
alignments with the major axes of field galaxies allowed us 
to include some of these in our sample. In MS 0451, we ob- 
served most galaxies with aligned slits, and, as a result, our 



field sample of spirals is weighted toward galaxies in the field 
around MS 0451. 

As our field sample is composed of objects in the same 
redshift range as the clusters, biases introduced by the 
magnitude-limited survey should affect both samples equally, 
except at the faint end of the luminosity function. This effect 
will be discussed further in § 4. 
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Fig. 3. — Continued. 



2.4. Source Extraction and Photometry 

Photometry was meas ured using SExtractor version 2.2.2 
dBertin & Arnou ts 1996). For ground-based imaging, we use 
SExtractor in two-image mode with source detection per- 
formed on the ground-based I-band images. Source detec- 
tion on the HST images was performed independently, and 
then matched to the ground-based catalog. For all imaging, 
we adopt magnitudes from the MAG_AUTO measurement of 
SExtractor. CI 0024 K$ images are not deep enough to detect 
all cluster spirals of interest; we therefore report K$ magni- 
tudes only for objects detected with > 3a significance. 

The kcorrect software v.4_l_2 (Bla nton et alj 120031) was 
used to convert observed K s and F814W fluxes to absolute 
magnitudes, Mk s (hereafter Mk) and My, respectively. At 
the redshifts of the clusters, observed /-band corresponds 
closely to rest-frame V, and so the k-corrections are small 
(< Q.5mag). In determining the k-corrections, we make use of 



all available ground-based photometry for each cluster field. 
This yields a better k-correction due to better sampling of the 
galaxy spectral energy distributions. For field objects at the 
low end of our redshift range, observed /?-band is a closer 
match to rest-frame V than observed /, and so for these ob- 
jects we apply a k-correction to the R magnitudes to deter- 
mine My- Before determining the k-corrections, we correct 
for a Galactic extinction of E(g-V)=0. 056 or E(g-V)=0.033 , 
for CI 0024 and MS 045 1 , respectively dSchlegel et al.|[l998b . 
All absolute magnitudes are expressed on the AB magnitude 
system. 

3. ROTATION CURVE ANALYSIS 

In order to construct the Tully Fisher relations, and to study 
other kinematic properties of spirals (such as mass, density, 
or M/L), we seek to determine the maximum velocity of the 
rotation curve for each disk. Our process involves extracting 
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the observed rotation curve from the spectrum and then creat- 
ing artificial rotation curves for each galaxy, determining the 
best-fit maximum velocity by x 2 fitting against the extracted 
rotation curve. In order to do this is important to determine 
various parameters about the galaxy from photometric data 
including: the position angle of the slit with respect to the po- 
sition angle of the major axis of the galaxy, the scale length 
of the galaxy and the seeing when the galaxy was observed. 
Along each step, we filter out the galaxies with weak spectral 
lines or those which, for reasons of inclination, etc. would 
not be possible to fit using our model. We further divide the 
remaining fits into two subsamples, those with secure rota- 
tion velocities (Q = 1), and those where the velocity is less 
certain but probably correct (Q = 2). We largely follow the 
procedure of Bohm et al. (2004), thoug h several other authors 
have followed similar procedures (e.g. iMetevier et all l2006t 
iNakamura et al.ll2006tlBamford et alj|2005UVopt et all 19961) . 
We have made several modifications to the procedure, detailed 
below. 

3.1. Extraction of Rotation Curves 

From each complete 2D spectrum, we extract postage 
stamps about the position of every emission line present, us- 
ing the known redshift of each object to identify lines. In the 
top panel of Figure Q] we display an example postage stamp 
centered on Ha for the galaxy N57426, a z = 0.39 field galaxy 
in the vicinity of MS 045 1 . As the observed center of Ha 
emission clearly varies across the spatial dimension of the 
spectrum, the rotation in this galaxy is already apparent. In 
order to determine the observed rotation curve, we fit a Gaus- 
sian function or double Gaussian (in the case of [Oil]) to each 
row along the spectral direction, as demonstrated in the lower 
panel of Figure Q] For the double Gaussian, we assumed a 
fixed separation and that the FWHM of each Gaussian com- 
ponent was the same, but allowed the amplitudes to differ in- 
dependently. 

We bin spectral rows together in the spatial dimension, as 
necessary, to meet a signal to noise (S/N) requirement of 

~ 5(A ) per bin. This re-binning allows us to sum up re- 
gions of the emission line that are too faint to fit individually, 
essentially trading spatial resolution, which is less important 
near the outer flat regions of a rotation curve, for a more re- 
liable velocity measurement. When rows are binned, the x 
position (i.e. radius) of the resulting velocity point is calcu- 
lated by taking a S/N weighted mean of the x positions of all 
the constituent spectral rows. 

Each of these binned fits were checked by hand to ensure 
that they were meaningful. All extracted rotation curves for 
our complete sample are displayed in Figures [2] and [3] 

We found that the fit worked best when we subtracted off 
the continuum to all the rows before fitting. We measured the 
continuum on the spectrum for each spectral line separately 
by summing together about 50 spectral columns on either side 
of the emission line's center. By extracting and continuum 
subtracting over small postage stamps, we avoid any issues 
that might be caused by spatial distortion in the spectrum. 

3.2. Surface Photometry 

From our HST imaging, we extract a 9" x 9" postage stamp 
image of each galaxy in our sample, with the galaxy centered 
See Figure s [2] and [3]). We then use the GALFIT software 
Peng et al.l2002l) to fit the galaxy photometry to a Sersic pro- 
file. This fit yields estimates of the galaxy position angle on 



the sky, axis ratio, and scale length. These parameters are 
necessary to determine the maximum velocity of the galaxy 
in the following step. While we experimented with fitting to 
an exponential disk plus bulge component, we find that the 
scale length as derived from the Sersic profile fit yields the 
best estimate of the rotation curve turn-over scale length, as 
described below. 

In local Tully-Fisher studies, internal extinction of spiral 
galaxies has been fo und to depend on the galaxy's inclination 
to the line of si ght (iTuilv & Fouqudll985t ITullv etal.|[l998l: 
IVerheiienl l2001). Galaxies viewed close to edge-on have a 
larger fraction of their luminosity extinguished by dust than 
the same galaxy would have if viewed face on. We correct 
for this effect by adoptin g the partic ularly simple form of the 
correction introduced by Tull v et all d!998l) : 

A A =- 7A log(^) (1) 

where a/ bis the axis ratio of the galaxy. This formula cor- 
rects toward the face-on case, but does not correct for addi- 
tional extinction in a face on galaxy. We choose not to apply 
any additional correction for the internal extinction of a face- 
on galaxy. 

Following rfullv et alj (1998), we determine 7 in each band 
by minimizing the scatter in the rest-frame color magnitude 
relations B-Ks vs. Mk and V -K$ vs. Mk, using the entire 
cluster plus field sample together. Since the luminosity func- 
tion of cluster spirals may not be uniform across all studied 
environments, we ignore any luminosity dependence of the 7 
correction, in order to avoid 'fitting away' real deviations that 
may be due to the cluster environment. We find -fs = 1.37, 
7v = 1.12, and jk = 0.15 x 7^ = 0.206. These values are con- 
sistent with the r ange specified in, e.g jTullv et al.l d 19981) and 
IVerheiienld200l . 

3.3. Model Fitting 

In order to determine the peak rotation velocity of a galaxy, 
we used the parameters obtained from GALFIT to construct 
an estimated velocity field for some maximum velocity. We 
adopt a standard rotation curve function of the form 

( r a + r ay/a 

where r s is the Sersic profile scale length a s dete rmined by 
GALFIT and a = 5, following iBohm et alJ (120041) . For two 
galaxies in our sample, pli97c4 and N46608, rs from GAL- 
FIT appeared to be an overestimate of the rotation curve scale 
length, and so we manually adjusted it to achieve a better fit. 
From the intrinsic rotation curve specified above, we construct 
a 2D velocity field by populating a grid of line-of-sight veloc- 
ities, under the following formula: 

V ohs = V(r)-cos(<f>)-sin(f> (3) 

where <j> is the azimuth in the plane of the disk and i is the 
inclination, with i = 90° defined to be edge-on to the line of 
sight. 

We then convolve this velocity field with a point spread 
function (PSF) with FWHM equal to the seeing. For our data, 
we used a fixed seeing of 0.8", equal to the median seeing of 
our observations. We adopt this fixed seeing correction be- 
cause of the relative insensitivity of the results to small varia- 
tions in seeing; we find that our uncertainty in the seeing cor- 
rection affects the final Vmax by < 1%, which is insignificant 
compared to errors due to inclination or position angle. 
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TABLE 1 
Summary of sample selection. 





CI 0024 


MS 0451 


Field 


Total 




0.37<z<0.41 


0.52 < z < 0.56 


0.3 <z<0.65 




a. Spirals in the specified redshift range with HST imaging and DEIMOS spectra 


103 


130 


194 


427 


b. Those with extended emission lines and aligned slits 


92 


103 


62 


257 


c. Significant spatial extent, with a measured velocity gradient 


42 


26 


47 


115 


d. After removing very face-on galaxies and other mis-aligned slits 


28 


16 


38 


82 


e. Velocity uncertainty small enough 


24 


16 


37 


77 


f. Q=l rotation curve 


17 


11 


30 


58 



Then, comparing the slit PA of our observation to the GAL- 
FIT estimate of the galaxy major axis angle, we place a mock 
l"-wide slit across the model velocity field, at an angle re- 
flecting the alignment between the real slit and the galaxy 
major axis. Finally, at each position along the length of the 
slit, we averaged the pixels across the slit width to determine 
an observed velocity. We use a % 2 minimization technique 
to vary the maximum velocity in the model to match our ob- 
served rotation curves. 

For each observed spectrum, we estimate the position of 
the galaxy's spatial center by fitting a Gaussian function to the 
spatial profile of the 2D absorption spectrum, integrated along 
the spectral dimension in two bands bracketing the emission 
line of interest. An initial estimate of the velocity center is 
calculated from the previously-determined redshift of each 
galaxy. Both the spatial center and velocity center are left 
to be free parameters in the \ 2 minimization. However, in no 
case does the best-fit spatial center of a Q = 1 rotation curve 
differ by more than two pixels (~ 1.5 kpc physical) from the 
calculated position, with typical offsets of much less than one 
pixel. 

We first fit the rotation curve using fixed values for the PA 
and inclination, i. In order to determine the error values on 
our fits, we factor in the error from the fit as well as comput- 
ing a PA error and inclination error for the model by running 
it at ±10° for each parameter. Especially for galaxies that 
present a somewhat face-on profile, it is important to account 
for this uncertainty due to PA and inclination errors, as it can 
be large in some cases, and in fact causes us to discard several 
emission line galaxies from our sample. We choose to vary 
over ±10° because the formal errors in the photometric fit 
are small in comparison to the systematic uncertainty in mea- 
suring the PA and inclination from the inherently asymmetric 
light profile of a spiral galaxy. 

3.4. Quality Control 

We began with a sample of 257 candidate spiral galaxies, 
each with visible, spatially resolved emission lines. Out of 
this sample, we removed 142 galaxies because the rotation 
curve lacked enough spatial extent to detect a reliable turn- 
over, or else no significant velocity gradient was measured, in 
most cases because the galaxy is oriented nearly face-on. Ad- 
ditionally, we removed 33 galaxies because the spectroscopic 
slit was too misaligned, or because the galaxy appeared too 
face-on to estimate the direction of its major axis. 

After fitting models to the rotation curves of our candidate 
spirals, and culling bad fits from our sample as described 
above, we remove five additional objects with highly uncer- 
tain velocities, ALog(2VMAx) > 0.2. Our final sample then 
consists of 37 field spirals and 40 cluster spirals (24 from 
CI 0024). The observed rotation curves, model fits, and im- 



ages of these galaxies are presented in postage stamp form in 
Figures [2] (cluster galaxies) and [3](field galaxies). 

We further divide this sample into two quality classes: Q=l 
rotation curves have turnovers detected on both sides of the 
curve, with a model fit that accurately matches the turnover at 
each end. Q = 2 curves, about 25% of the total, only show a 
turnover at one end, or show other signs of an uncertain fit to 
the model. We do not simply throw out all galaxies with signs 
of disturbed kinematics, but rather keep them in the Q = 2 
sample, as they are of considerable interest for our study of 
possible interactions with the cluster environment. However, 
the requirement that we identify a reasonably secure value of 
Vmax must necessarily exclude some number of spirals with 
highly disturbed rotation curves. 

We summarize the sample selection process for both clus- 
ters and the field in Table 1 . We note that we remove roughly 
equal fractions of galaxies from the cluster and field samples, 
at each step of the process. The two notable exceptions have 
ready explanations: In step b.), a larger fraction of field galax- 
ies than cluster galaxies are removed due to misaligned slits, 
because we did not consistently observe field galaxies with 
aligned slits. Similarly, in step c), a low fraction of MS 0451 
cluster galaxies exhibited large enough spatial extent in their 
rotation curves, due to an observed suppression or lack of star 
formation across this massive cluster (Moran et al. 2007, in 
prep). Basic data, as well as extinction-corrected magnitudes 
and velocities for all 77 objects of our main sample are listed 
in Table 3. 

4. THE TULLY FISHER RELATION 

In Figure |U we plot the Tully Fisher relation for both clus- 
ter and field galaxies, in both rest frame K$ and V bands (ex- 
pressed as absolute magnitudes Mk and My, respectively). 
Shad ed regions indica te the ± 1 - a scatter of the local TF 
from IVerheijenl (120011) . Solid lines indicate the best fit TF 
zero point for each o f our subsamples , with the slope fixed to 
the local values from Verheijen] d2001l) . adopting their RC/FD 
sample, which includes only galaxies where a turnover in the 
rotation curve is seen. Dotted lines indicate the 1 — a RMS 
scatter of the relation about the mean. Zero points are cal- 
culated by finding the bi-weight mean of the residuals about 
the local TF. To minimize bias, zero point and scatter are cal- 
culated only from Q= 1 rotation curves. We further impose 
an absolute magnitude cut of My < -20.2 and Mk < -20.2, 
to eliminate bias due to the differing magnitude distributions 
between the cluster and field samples. 

For purposes of determining the intrinsic scatter in the TF 
relation for each of our samples, we also perform a least- 
squares fit to find the best-fit parameters of the TF rela- 
tion, weighting each point by the measurement uncertain- 
ties in both Log(2VMAx) and absolute magnitude. We follow 
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FIG. 4. — V-band (right) and K-band (left) Tully-Fisher relations for cluster and field spirals. The field relation is plotted in the top two panels, and cluster 
galaxies are plotted on the bottom. CI 0024 members are marked as black circles, while MS 0451 members are indicated by red diamonds. Shaded regions 
indicate the 1 — a scatter of the local Tully-Fisher relation, from Verheijen (2001). We adopt their R-band relation to compare to our V-band relation, neglecting 
any small correction to account for the different filters. In each panel, the thick solid line represents the inverse-fitted TF relation for those points (but fixing 
the slope to the local value), and the dotted lines represent the la scatter about the mean. Dash-dotted line indicates where we apply a cut in magnitude for this 
calculation. Symbols are as indicated by the legend; small symbols indicate Q = 2 rotation curves, and large symbols are Q = 1 . 



other authors (e .g. Bamf ord et aT1l2006t iMetevier et al.ll2006t 
iNakamura et al.l 120061) and adopt Log(2VMAx) as the depen- 



dent variable in the fit, such that 



Log(2V MAX ) = a+b*M v 



(4) 



where we fit for intercept a and slope b. This is the so-called 
'inverse TF', w hich is less sensitive to bias du e to luminosity 
incompleteness ( Willickll 1 994t ISchechterll 19801) . For each TF 
relation, cluster and field, V and K s bands, the intrinsic scatter 
is that portion of the measured scatter that cannot be explained 
by measurement error. We estimate the intri nsic scatter by 
considering the reduced \~ statistic. Following Ba mford et al.l 



(2006), we iteratively determine the scatter that we need to 
add to our measurement errors in order to achieve x? = 1 • The 
best inverse fit parameters for all four subsamples are listed 
in Table 2. We note that the zero points and slopes are in- 
distinguishable between cluster and field; as we will discuss 
below, only the scatter in the relation differs between cluster 
and field. 

Field galaxies show a tight TF relatio n in both bands, with 
slope consistent with the local relation (Verheiien 200l]). We 
find an intrinsic scatter of 0.35 mag in V and 0.5 mag in Ks, 
again restricting ourselves to Q = 1 rotation curves brighter 
than our magnitude cut. The seemingly higher intrinsic scat- 
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TABLE 2 

Inverse fits to Tully-Fisher relation. 



Sample 


a 

RMS 


b 

RMS 


(ALog(2V MAX )) 


(AM K . V ) 


Field K s 


-0.7 ±0.3 


-0.14 ±0.01 


0.07 ±0.01 


0.5 ±0.05 


Cluster K s 


-0.5 ±0.2 


-0.14 ±0.01 


0.17 ±0.01 


1.18±0.05 


Field V 


-2.1 ±0.4 


-0.22 ±0.02 


0.08 ±0.01 


0.35 ±0.05 


Cluster V 


-1.9±0.2 


-0.20 ±0.01 


0.19±0.01 


0.93 ±0.05 



ter in K s is at odds with the expectation that lower dust extinc- 
tion in the Ks-band should make the TF tighter than in bluer 
bands, and our result seems to indicate that we underesti- 
mate the measurement uncertainties. However, the shallower 
slope of the Ks-band TF causes its scatter, when expressed in 
magnitudes, to be more sensitive to small errors in the mea- 
sured Vmax- In fact, expressed in terms of ALog(VMAx), the 
field sample scatter in V and Ks are indistinguishable, yield- 
ing 0.08 ±0.01 and 0.07 ±0.01, respectively. As measure- 
ment uncertainties are generally largest in the Vmax direction, 
this simply indicates that absolute magnitude is more prop- 
erly the independent variable in our TF relation. In the fol- 
lowing we will preferentially express the TF scatter in terms 
of ALog(2VMAx), except when comparing to other authors. 

The 0.35 mag scatter we find i n V is comparable to the 0.38 
mag R-band scatter reported by IVerheiie n (2001) for nearby 
galaxies. Our A^-band scatter of 0.5 mag is about 50% larger 
than t heir reported 0.31 mag. However, iKannappan et alJ 
(2002) have suggested that, because local studies tend to weed 
out kinematically irregular galaxies, the true scatter, if a more 
representative sample of spirals is selected, could actually be 
much higher. At intermediate redshift, small irregularities in 
rotation curves are harder to detect, due to limited spatial res- 
olution, and so a higher measured scatter might reasonably be 
expected. 

Even so, the scatter we measure for the field TF is signifi- 
cantly lower than h as been found by other authors at interme- 
diate redshift (e.g. iNakamura et ai1l2006t iBohm et alJ 120041: 
IBamford et al.ll2.Q06h . In a recent large study of 89 field galax- 
ies, iBamford et alJ ([2006) measure an intrinsic scatter in the 
B-band TF of 0.9 mag, significantly larger than our measured 
V-band scatter, though they include galaxies across a larger 
redshift range. In the redder bands that we measure, we can 
see that a tight TF relation still exists in the field at lookback 
times of over 5 Gyr. 

In stark contrast to the field TF, the relation that we measure 
for cluster spirals (lower panels of Figure |4| shows a remark- 
ably high scatter, ALog(V M Ax) = 0.19 ±0.01 (0.93 mag) and 
0.17±0.01 (1.18 mag) in V and Ks, respectively. These val- 
ues are each more than twice as large as the scatter in our 
field sample. This cannot be understood in terms of higher 
measurement error in the cluster sample, as the two samples 
were selected in the same way from the same parent data set, 
and we have restricted the analysis to only the highest qual- 
ity rotation curves. When we include all 77 rotation curves, 
we observe the same difference between cluster and field, but 
with overall higher measurements of scatter. 

4.1. A Comparison to Independent Measurements 

Recently, iMetevier et alJ d2006) have also published a 
Tully-Fisher relation for the cluster CI 0024, examining rota- 
tion curves of 15 spirals. Four of the galaxies in their sample 
are in common with our own, allowing us for the first time 
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FIG. 5. — The residuals from the local V-band Tully-Fisher relation, in 
the log(2V M/ tx) direction, versus rest (V -Ks)ab color. Magnitudes are cor- 
rected for inclination-dependent internal extinction, as discussed in the text. 
Symbols are coded as in Figuref4] 

to evaluate the agreement between repeat observations and 
independent analysis of intermediate redshift spiral rotation 
curves. In our sample, galaxies p0i27c3, p0il45c4, p0il63c4, 
and pl9i2c4 correspond to their TFR05, TFR07, TFR10, and 
TFR12, respectively. 

Visual comparison of their observed rotation curves to our 
own (Figure [2]) indicates that they are of comparable quality, 
but with some differences in rotation curve extent. Comparing 
our estimates of Vmax to their Varc, we find an RMS differ- 
ence of AV jV = 34% ± 17%, with individual measurements 
differing by as much as 100kms~ l . Three out of four measure- 
ments differ by more than 2cr. 

In IMetevier et al. (2006), Varc is conceptually identical to 
our Vmax- both attempt to measure the broad flat part of each 
spiral's rotation curve. Furthermore, our procedure for mod- 
eling the rotation curve follows steps very similar to their 
GAUSS2D code. However, we each adopt slightly differ- 
ent rotation curve functions; they use an arctan function to 
approximate V(R), while we use the fu nction given in Equa- 
tion 2, adopted from Boh m et alJ (120041) . 

To test the effect of adopting a different rotation curve func- 
tion, we re-run our model fi ts for the four galaxies in common 
with [Metevier et al. (2006), this time fitting the observed ro- 
tation curves to an arctan function, and adopting the best-fit 
scale lengths from their paper. We find that adopting their ro- 
tation curve function brings two of our four velocity measure- 
ments into agreement. For the two other objects, variations in 
the observed rot ation curves may ex plain the discrepancy. In 
p0il63c4/TF10, Met evier et al.1 d2006l) uncover a downturn in 
the rotation curve at high radius, which is not reached by our 
own data. Conversely, for object p0il45c4/TF12, our rotation 
curve extends to larger radius and reveals that the velocity 
co ntinues to increase be yond the end of the curve measured 
bv IMetevier et all d2006l) . 

It is striking that such a small difference in the choice of 
rotation curve function can yield such a large difference in 
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FIG. 6. — The residuals from the local Tully-Fisher relation as a function of normalized projected radius, R/Rvir. Symbols are coded as in Figure|4] Black 
crosses to the right of the dash-dotted line indicate the median and RMS residuals of the field TF. Positive residuals indicate that a galaxy is overluminous for its 
measured velocity, or, alternatively, has an anomalously low Vmax given its luminosity. 



the resulting velocity value . When compar i ng the ir results to 
other Tully-Fisher studies, iMetevier et alj (|2006) take pains 
to comprehensively account for many of the systematic dif- 
ferences between studies, most having to do with the way 
velocities were measured or defined. They show that these 
small differences greatly affect estimates of the average lumi- 
nosity evolution of spirals as a function of redshift. Because 
of the difficulty in comparing TF relations across samples, 
and the additional, previously under-appreciated systematic 
arising from the choice of rotation curve function, we do not 
attempt in this paper to make a rigorous estimate of the lumi- 
nosity evolution implied by our TF zero point. 

In fact, generally, and for our present goal of studying en- 
vironmental influences on spiral galaxies, these variations be- 
tween studies highlight the importance of matched samples of 
cluster and field spirals. Our current matched sample, which 
contains the largest number of cluster galaxies so far, allows 
us to move beyond the Tully Fisher relation, to study directly 
the kinematics of spirals galaxies in the hope of uncovering 
the source of the large scatter seen in the cluster TF relation. 

5. TRENDS IN SPIRAL MASSES, DENSITIES, AND M/L 

Movement of a galaxy in the TF plane can be caused by 
various effects. Increased star formation or dustiness moves a 
galaxy along the luminosity axis, while increased total mass, 
changes to the radial mass profile, or other kinematic distur- 
bances can alter a galaxy's measured Vmax- In this section, 
we attempt to identify the source of the high scatter observed 
in the cluster TF, compared to the field in the same redshift 
range. Because our cluster sample spans a large range in 
cluster-centric radius, it makes sense to examine the residuals 
of the TF and other dynamical characteristics of the cluster 
spirals as a function of environment within the clusters. 

5.1. Tully Fisher Residuals and M/L 

We first turn toward the simplest quantity to examine: the 
residuals from the TF, considered to represent a change in 



mass to light ratio (M/L). The higher scatter seen in the cluster 
TF, compared to the field, could plausibly be due to changes 
in star formation rate or dust obscuration as galaxies fall into 
the cluster. This could cause a radial gradient in M/L, which 
we would observe as an increased TF scatter. 

However, because enhanced scatters are measured in both V 
and K$ band TF relations, we do not expect that the enhance- 
ment can be solely attributed to an increase or decrease in 
star formation or dust during cluster infall. If such a scenario 
were the main driver of the scatter, we would expect that the 
V-band scatter would be higher than the /^-band, yet they are 
broadly equivalent. We would also expect to see a correlation 
between TF residuals and other indicators of star formation 
rate or dustiness. To test this, we plot, in Figure[5] the V-band 
TF residuals versus V— Kg color for all galaxies in our cluster 
sample. Here, and for the rest of this paper, residuals are plot- 
ted in the sense that a positive residual is overluminous for 
its measured velocity, or, alternatively, has an anomalously 
low Vmax given its luminosity. No obvious correlation is ob- 
served in Figure [5] While several of our objects with lowest 
measured Vmax are not plotted here because they do not have 
Ks detections, these were also excluded from the TF relation 
fitting, and so are not the source of the enhanced scatter. 

Even though star formation rate and dust content do not ap- 
pear to correlate with TF residuals, the galaxies of our cluster 
sample span a wide range of environments, and so we also ex- 
amine whether the increased TF scatter is related to a gradient 
in M/L across the cluster, as might be the case if spirals in the 
outskirts have formed more recently than those in the clus- 
ter cores. In Figure [6] we plot the TF residuals as a function 
of R/Rvir, the projected cluster-centric radi us scaled by each 
cluster's Virial radius: 1.7 Mpc for CI 0024 dTreu et al.ll2003l) 
and 2.7 Mpc for MS 0451 dMoran et al.ll2007l) As MS 0451 
is much more massive than CI 0024, a galaxy at 1 Mpc radius 
in CI 0024 experiences a very different environment from a 
galaxy at the same radius in MS 045 1 . As we will see below, 
some key trends emerge when we choose to scale by Virial 
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FIG. 7. — Left: Projected central density of cluster spirals, measured within 1.25r s , versus projected R/Rvir, calculated according to £,„(< r) = V 2 /(Gr) Right: 
Central density versus radius R. In both panels, The thick cross to the right of the dash-dotted line indicates the median density for our field sample. In the left 
panel, additional thick crosses indicate the median densities inside and outside 1 Rvir- Symbols are coded as in Figure|4] 

radius, rather than plotting radius directly on the x-axis. 

Examined by eye, the residuals from the V band TF in Fig- 
ure [6] (right panel) hint at a radial gradient as a function of 
R/Rvir, with galaxies at higher radius seeming to be over- 
luminous. However, straight line fits to the V and Ks band 
residuals find a small gradient toward higher radius, but with 
slope no greater than the error bar on a typical point. A sim- 
ple gradient in star formation rate or M/L across the cluster 
therefore cannot be the only mechanism responsible for the 
increased scatter in the cluster TF compared to the field. We 
note, however, that we cannot rule out the possibility that sev- 
eral different mechanisms are simultaneously contributing to 
the TF scatter by acting on spirals in different environments 
within the clusters. 

5.2. Densities and Masses 

Since variations in star formation rate and dust content 
alone cannot account for the observed scatter in the cluster TF, 
we are led to consider the idea that the cluster spirals are more 
kinematically disturbed than their field counterparts. One way 
to test for disturbed dynamics in a spiral galaxy is to consider 
the photometric effective radius, r s , of each galaxy. We can 
combine r s with Vmax to calculate two fundamental dynami- 
cal properties of spiral disks: dynamical mass, M(< r) oc V 2 r 
and central surface mass density, S m (< r) ocV 2 /r. Unlike the 
Fundamental Plane of ellipticals, the Tully-Fisher relation in 
the local universe d oes not seem to have any dependence on 
galaxy size (r s ) (e.g. lVerh eiien 2001). Therefore, in an undis- 
turbed population of spirals, we would not expect to uncover 
any independent environmental trends in quantities that only 
depend on r s and Vmax- Any trends that do exist must be the 
result of some cluster-related physical process. 

In fact, surface densities allow us to directly probe for 
the action of a key physical mechanism, galaxy harassment 
dMoore et al.l Il999h . Harassment is predicted to have a 
stronger effect on the least dense galaxies falling into a clus- 
ter, to the point of completely disrupting the most tenuous 
spirals. Any observed gradient in the mean density of spirals 




FIG. 8. — Dynamical mass (within 1.25r,) as a function of R/Rvir, cal- 
culated according to M(< r) = V 2 r/G. We find galaxies with a range of 
dynamical masses, both within and outside of the Virial radius. Symbols are 
coded as in Figure[4] The median for the field sample is displayed to the right 
of the dash-dotted line. 



could then implicate the action of this physical process. 

We choose to study M and £,„ within a radius of 1.25r s , a 
characteristic radius chosen because it is typically reached in 
all of our observed rotation curves, and it is a radius at which 
most of our rotation curves have already leveled off to Vmax- 

In Figure [7] we plot galaxy density as a function of pro- 
jected radius R (right panel) as well as the normalized quantity 
R/Rvir (left). In the left panel, one notices a striking break in 
the densities of spirals at approximately IRvir- Near and out- 
side of this radius, spirals seem to exhibit nearly uniformly 
low central densities, which are puzzlingly even lower than 
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those of field galaxies in the sample. Within the Virial ra- 
dius, on the other hand, a large spread in densities is seen, 
and perhaps a radial gradient of decreasing density outward 
from the cluster center. In the right hand panel, this trend 
appears scrambled, indicating that whatever physical process 
may cause this effect, its strength scales as the cluster Viral ra- 
dius. This observation rules out several possible mechanisms, 
and will be discussed further in the next section. 

It is natural to wonder if the observed break in density as 
a function of radius is due to simple luminosity segregation: 
if more massive spirals are found near the cluster center, per- 
haps these also have higher central densities. However, by 
consulting Figure[S] where we plot dynamical mass as a func- 
tion of R/Rvir, it becomes apparent that this is not the case. 
The low density spirals found in the clusters' outskirts in fact 
exhibit a wide spread in total mass. Because the total number 
of observed spirals in the cluster outskirts is small, K-S tests 
comparing the distributions of mass and density for low-R vs. 
high-R spirals are inconclusive. However, the cluster sam- 
ple as a whole does exhibit a larger overall spread in £,„ than 
the field sample, and includes a larger fraction of both low 
and high density spirals: a K-S test gives a 1.3% chance that 
the two samples are drawn from the same parent distribution. 
This indicates that the cluster environment may be affecting 
the internal mass distributions of spirals at all cluster radii. 

6. DISCUSSION 

What physical mechanisms, then, could be acting on cluster 
spirals to reproduce both the overall higher scatter in the TF 
and the observed radial trend in density? The effects seem to 
persist as far as 2Rvir from the cluster cores, so even though 
some of the observed galaxies at high radius may be part of 
a "backsplash population", it is very unlikely that nearly all 
star forming spirals in our sample have already been through 
the cluster center. Therefore, we do not think it likely that 
tidal processes are responsible, as they are only strong near 
the cluster center. Instead, we can consider several proposed 
physical mechanisms that are strong enough at large radius to 
alter the dynamics of a spiral, either directly or indirectly. 

In the cluster outskirts, recent mergers are an obvious candi- 
date to drive bo th the large TF scatter and the abnormally low 
densities (e.g. lTreu et al.| [2003, and references therein) Merg- 
ers are less common in the low density field, and this could 
explain why we do not see any low-density spirals in our field 
sample. Unknown selection biases, though, could prevent us 
from including similar field galaxies in our Tully-Fisher sam- 
ple. If the ef fects of a recent merger last for at least 1 Gyr 
(Bekkii [T998l) . then rece nt mergers can affect the TF scatter 
even in the cluster core dTreu et al.ll2.Q03b . Therefore, we can 
not rule out the possibility that increased merging in the clus- 
ter outskirts serves to drive a high fraction of cluster spirals 
away from the TF relation. 

In the inner regions of galaxy clusters, mergers are sup- 
pressed due to the high relative speeds of galaxies, which pre- 
vent the creation of a gravitationally bound pair during close 
encounters. Instead, an infalling cluster galaxy is likely to ex- 
perience repeated close encounters at high speed due to the 
high density of galaxies in the cluster. This process, called 
galaxy- galaxy haras s ment, can lead to dramatic changes in a 
galaxy. iMoore et al.l i 1999b have shown through simulations 
that the fate of a harassed galaxy depends on its original mass 
and central density. Strongly concentrated Sa/Sb type galax- 
ies were seen to puff up their disks during infall, and so ha- 
rassment may represent one way in which spirals transform 



into SOs in clusters. On the other hand, lower d ensity Sc/Sd 
spirals are more strongly affected by harassment lMoore et al.1 
(1999) found that they' were either completely disrupted, or 
else transformed into an object resembling a dwarf galaxy. 

If harassment is acting to transform the lowest density spi- 
rals into dwarfs, then we would expect to observe a deficit of 
such low-density spirals near the cluster cores. High-density 
spirals, on the other hand, should be more resistant to harass- 
ment, and are likely to persist to smaller cluster radii. This 
prediction qualitatively matches the trend in densities seen in 
Figure|7] but the picture is unclear. The puzzling lack of high- 
density galaxies at large cluster radius, and the persistence of 
low-density galaxies to ~ Q.SRvir raise questions about this 
interpretation. We have already seen in Figure [8] that galax- 
ies of a wide range of masses are represented in the cluster 
outskirts, so harassment alone may not present a complete ex- 
planation for the observations. Frequent mergers in the cluster 
outskirts, however, could very well provide the missing ingre- 
dient for keeping spiral densities low there. 

Finally, we consider processes that depend on the hot in- 
tracluster medium (ICM). Generally, even strong interac- 
tions with the ICM like ram pressure stripping are thought 
to be too weak to explain the observed di sruptions in the 
kinematics of spiral disks (lOuilis et al.l l2000h . Rather, such 
ICM-related processes act largely to suppress star formation 
within infalling disks. Since changes in star formation rate 
do not appear to be responsible for the increased scatter in 
the cluster TF, it is unlikely that an ICM-related process is 
involved. However, it is possible that shock fronts within 
the ICM can enhance the IC M's ability to affe ct a spiral 
disk even at high cluster radius dMoran et al.ll2005l and refer- 
ences therein). CI 024 may be undergo ing a face-on merger 
with a large group (ICzoske et al.1 12002). and so shocks may 
be important in this cluster. Shocks in the ICM may in- 
duce cent rally-concentrated starbursts within infalling cluster 
galaxies dMoran et al.ll2005l) . but it is not clear that such an 
interaction would generate emission lines with enough spa- 
tial extent to allow measurement of rotation curves. Further, 
since all ICM related processes suppress star formation over 
time, our sample of exclusively star-forming spirals can only 
provide and incomplete picture at best of the effects of ICM 
shocks. 

One possible concern with our result on the cluster TF re- 
lation is that cluster to cluster variation could be high (as 
seen for example by the MORPHS a nd EDIScS studies, 
(Pogg ianfi et"atlll999l: IWhite et al.ll2005l) ). In fact, previous 
studies of CI 0024 have shown that its galaxy population 
may be overly active, possibl y due to the ongoi ng merger 
with a foreg r ound l arge group dCzoske et al.ll2002h . Indeed, 
iMoran et al.l (120071) showed that the Fundamental Plane of 
elliptical and SO galaxies in CI 0024 exhibits a higher scat- 
ter than found in m ost other intermediate redshift clusters 
dKelson et al.ll20061 e.g.), though this effect is most signifi- 
cant in the inner 1 Mpc of the cluster. It is possible that the 
increased TF scatter we see is connected to the similarly en- 
hanced Fundamental Plane scatter. However, while our clus- 
ter TF includes a majority of points from CI 0024 (24 galaxies 
c.f. 16 for MS 0451), it seems clear by inspection of Figure[4] 
that MS 0451 also contains spirals that deviate highly from 
the local TF. As MS 045 1 is thought to be in a more advanced 
stage of cluster assembly than CI 0024 (Moran et al. 2007, in 
preparation), the universality of our measured TF scatter re- 
mains uncertain until similarly large samples for several more 
intermediate redshift clusters become available. 
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7. CONCLUSIONS 

In this paper, we have studied the dynamics of cluster and 
field spirals at intermediate redshifts, via an analysis of their 
optical rotation curves. 

We have presented one of the most complete Tully Fisher 
relations available for both cluster and field spirals at these 
redshifts, and have demonstrated that the field relation is quite 
tight even at these high look back times. In contrast, the clus- 
ter TF exhibits a remarkably high scatter. By comparing the 
trends of Tully-Fisher residuals vs radius with those in col- 
ors and local mass density we found that the increased scatter 
cannot be explained solely in terms of environmental effects 
on the star formation rates of infalling galaxies. We therefore 
proposed that the increased scatter in the Tully Fisher relation 
is due to kinematic disturbances, as expected for example for 
cluster harassment. 

We also found a trend in galaxy mass density as a func- 
tion of cluster centric radius in the sense that spiral galaxies 
are denser in the central regions of clusters. This is expected if 
harassment plays a significant role, as the densest disks would 
be most likely to survive during the infall. However, we found 
a paucity of high density spiral galaxies in the cluster out- 
skirts, which cannot be explained by harassment alone. We 
suggest that a combination of enhanced merging in the cluster 
outskirts with galaxy harassment in the intermediate and inner 



cluster regions may be required to explain the observed trend 
in galaxy density. 

Larger matched samples covering a larger number of galaxy 
clusters are needed to determine if the observed trends are 
universal across clusters at these redshifts. 
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TABLE 3 

Information and Measurements of Cl 0024, MS 0451 and field galaxies 



Object 


Sample 


z 


RA 


DEC 


r s 


/ 


SI PA 


M v 


M K 


%AX Q 








[deg] 


[deg] 


[kpc] 


[deg] 


[deg] 


[mag] 


[mag] 


[km/s] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 



p()i27c3 


C10024 


0.392 


6.639747 


P 38il4c3 


C10024 


0.400 


6.631994 


p0il63c4 


C10024 


0.394 


6.674678 


p0il45c4 


C10024 


0.399 


6.676911 


N36671 


MS0451 


0.527 


73.521355 


N37826 


MS0451 


0.530 


73.523247 


p38il22c2 


C10024 


0.390 


6.618957 


N35977 


MS0451 


0.531 


73.506584 


P 37ill0c3 


C10024 


0.377 


6.671664 


P51837 


MS0451 


0.526 


73.566147 


p37il8c3 


C10024 


0.392 


6.676580 


P 36i29c3 


C10024 


0.400 


6.616654 


N25094 


MS0451 


0.535 


73.499481 


P62230 


MS0451 


0.530 


73.562500 


p36i77c4 


C10024 


0.396 


6.640263 


P63232 


MS0451 


0.538 


73.494949 


p37i2c3 


C10024 


0.377 


6.676347 


p36i37c4 


C10024 


0.388 


6.630331 


P54663 


MS0451 


0.541 


73.597534 


pl9i2c4 


C10024 


0.396 


6.591903 


pl9ilc4 


C10024 


0.397 


6.595224 


P 13il9c4 


C10024 


0.381 


6.617717 


P 23il51c2 


C10024 


0.393 


6.679020 


P 18i234c4 


C10024 


0.396 


6.721228 


N27099 


MS0451 


0.550 


73.635056 


N49410 


MS0451 


0.528 


73.533180 


P 24i55c2 


C10024 


0.388 


6.639475 


P 25il36c3 


C10024 


0.395 


6.584536 


p24ilc2 


C10024 


0.392 


6.646124 


N46608 


MS0451 


0.540 


73.473412 


N48527 


MS0451 


0.544 


73.488052 


pl3il5c2 


C10024 


0.397 


6.598421 


N34128 


MS0451 


0.534 


73.427948 


N48766 


MS0451 


0.551 


73.469215 


N47732 


MS0451 


0.534 


73.413109 


N64729 


MS0451 


0.540 


73.520599 


p7i72c2 


C10024 


0.373 


6.597922 


p30i4c3 


C10024 


0.393 


6.552849 


pl0i41c4 


C10024 


0.396 


6.801073 


p33i3()8c4 


C10024 


0.382 


6.603914 


P 5il71c4 


field 


0.313 


6.742421 


N33315 


field 


0.314 


73.373642 


N22367 


field 


0.325 


73.610954 


N47194 


field 


0.326 


73.525711 


N58982 


field 


0.326 


73.631409 


N11649 


field 


0.329 


73.513878 


N30997 


field 


0.333 


73.444008 


N11932 


field 


0.362 


73.463303 


N48819 


field 


0.363 


73.460930 


N55288 


field 


0.364 


73.419220 


N37713 


field 


0.367 


73.467995 


N21079 


field 


0.371 


73.350967 


N31950 


field 


0.371 


73.444237 


N35309 


field 


0.390 


73.625328 


N57426 


field 


0.391 


73.473892 


N18456 


field 


0.401 


73.598076 


N35519 


field 


0.413 


73.426918 


N 19992 


field 


0.419 


73.515671 


N18483 


field 


0.425 


73.465881 


P 17il29c4 


field 


0.443 


6.766013 


N36846 


field 


0.447 


73.670662 


N40959 


field 


0.447 


73.624619 


N41465 


field 


0.447 


73.629173 


N56286 


field 


0.463 


73.505577 


pllillc4 


field 


0.476 


6.724503 


N25636 


field 


0.491 


73.394081 


p8il63c3 


field 


0.492 


6.569071 


N59564 


field 


0.494 


73.547920 


pli97c4 


field 


0.494 


6.770693 


N11338 


field 


0.505 


73.540504 


N25124 


field 


0.506 


73.698349 


p23il80c2 


field 


0.536 


6.677425 


P 21i99c3 


field 


0.537 


6.475233 



17.156639 
17.154470 
17.164650 
17.158810 
-2.994348 
-2.988543 
17.155710 
-2.993443 
17.129601 
-3.065880 
17.128189 
17.189550 
-3.053076 
-3.073720 
17.205460 
-2.982716 
17.117929 
17.208891 
-3.063867 
17.168131 
17.186880 
17.218800 
17.103390 
17.167561 
-3.040075 
-2.924410 
17.082970 
17.114901 
17.081320 
-2.938749 
-2.927898 
17.232731 
-3.003162 
-2.925409 
-2.930732 
-2.835176 
17.299730 
17.050751 
17.199739 
17.000919 
17.238001 
-3.008564 
-3.070553 
-2.939083 
-2.867468 
-3.123671 
-3.020530 
-3.124372 
-2.927250 
-2.892889 
-2.987958 
-3.074297 
-3.017617 
-2.996757 
-2.878548 
-3.087556 
-2.997052 
-3.081446 
-3.088012 
17.151690 
-2.991703 
-2.969724 
-2.963893 
-2.884414 
17.193130 
-3.050711 
17.296320 
-2.870288 
17.324760 
-3.126556 
-3.052146 
17.101549 
17.202730 



5.23 
4.91 
8.70 
3.05 
6.68 
6.18 
2.38 
3.48 
3.82 
7.33 
6.50 
2.66 
6.24 
3.53 
7.91 
4.23 
10.23 
4.69 
3.35 
3.71 
3.94 
4.26 
3.02 
2.55 
4.67 
8.54 
1.85 
5.37 
3.19 
1.75* 
9.94 
9.00 
4.28 
3.49 
4.23 
4.57 
1.68 
5.82 
13.68 
7.10 
3.43 
4.16 
5.06 
7.08 
4.26 
2.40 
2.59 
3.95 
1.93 
5.20 
4.95 
2.86 
3.69 
3.32 
3.73 
4.13 
3.83 
6.52 
4.56 
10.24 
2.73 
3.40 
4.67 
4.59 
2.77 
6.04 
3.14 
4.66 
3.00* 
5.12 
7.84 
5.56 
2.68 



54.2 
50.5 
62.0 
28.2 
39.6 
57.4 
36.1 
67.4 
80.4 
62.3 
69.8 
79.8 
68.7 
55.2 
77.6 
47.4 
75.5 
69.2 
56.1 
31.0 
32.8 
78.8 
23.8 
62.9 
77.9 
35.4 
71.6 
81.9 
60.4 
61.1 
50.3 
75.9 
75.5 
33.9 
46.3 
44.2 
46.0 
35.8 
41.4 
46.8 
74.8 
55.8 
74.7 
33.0 
44.3 
57.4 
62.4 
41.0 
28.1 
62.6 
27.2 
42.9 
55.2 
66.3 
60.3 
61.0 
40.8 
51.5 
71.7 
84.6 
45.3 
34.8 
74.7 
53.6 
27.3 
28.5 
56.7 
71.4 
38.8 
36.2 
39.4 
56.0 
36.1 



11 

60 
3 
13 
36 
10 
41 
13 
22 
55 


29 
42 
44 

28 

1 

21 
29 
8 

21 
16 
29 
53 
27 
24 
13 
21 
40 

3 

26 
24 
38 
18 
54 
30 
26 

47 
45 
44 
14 
42 
20 
67 
26 
50 


33 
41 
21 
38 
14 
43 
4 
3 
2 

13 
21 
3 
18 
15 
7 
22 
16 
33 
44 
14 
3 



21 


80±0.05 


-22 


58±0 


06 


219 ± 14 




21 


67±0.06 


-22 


32±0 


06 


322 ± 119 


j 


22 


89±0.()9 


-23 


55 ±0 


05 


300 ± 10 


j 


20 


43±0.07 


-20 


05 ±0 


30 


147 ± 19 


j 


22 


89±0.07 


-23 


47±0 


05 


105 ± 16 




23 


12±0.06 


-23 


75±0 


05 


316 ± 18 


1 


19 


84±0.06 


-19 


84±0 


32 


61 ± 11 


j 


20 


30±0.07 








122 ± 7 


I 


20 


81 ±0.07 








78 ± 9 




21 


77±0.08 


-22 


04±0 


05 


317 ± 90 


j 


20 


75±0.05 








109 ± 3 


j 


20 


62±0.07 








45 ± 6 


j 


21 


56±0.07 


-21 


61±0 


06 


1 12 ± 20 




21 


21±0.06 


-21 


72±0 


05 


89 ± 16 


j 


22 


22±0.06 


-22 


02±0 


10 


143 ± 4 


j 


20 


92±0.06 


-21 


35 ±0 


06 


158 ± 21 


I 


21 


31±0.13 


-21 


76±0 


20 


167 ± 4 


2 


21 


07±0.06 


-20 


65 ±0 


23 


105 ± 2 


2 


21 


65±0.06 


-22 


35±0 


05 


152 ± 14 


2 


21 


47±0.06 


-22 


62±0 


05 


152 ± 26 




20 


69±0.07 


-20 


8()±() 


20 


86 ± 10 




21 


34±0.05 








96 ± 1 


j 


20 


32±0.05 


-22 


09±0 


06 


93 ± 16 




20 


14±0.08 


-20 


25±0 


29 


121 ± 6 


j 


21 


73±0.06 


-21 


35±() 


07 


189 ± 21 




21 


70±0.05 


-21 


93±0 


05 


150 ± 47 


1 


19 


29±0.09 








30 ± 11 


j 


21 


16±0.08 








67 ± 5 




19 


93±0.10 


-20 


85±0 


24 


26 ±4 




21 


79±0.05 


-22 


67±0 


05 


139 ± 13 


j 


22 


05 ±0.05 


-22 


60±0 


05 


341 ± 62 


I 


21 


29±0.05 








136 ± 4 




20 


84±0.06 








55 ± 9 


j 


21 


30±0.05 


-22 


14±0 


05 


87 ± 11 


j 


20 


26±0.06 


-20 


17±0 


15 


94 ± 11 


j 


21 


37±0.06 


-21 


8()±0 


05 


84 ± 12 




19 


04±0.07 








27 ± 3 


j 


21 


46±0.13 


-22 


48±0 


10 


398 ± 138 


j 


22 


30±0.10 


-23 


()4±0 


05 


110 ± 18 


j 


20 


21±0.05 


-20 


65 ±0 


25 


110 ± 13 


I 


20 


77±0.11 


-21 


38±() 


07 


99 ± 3 


j 


21 


67±0.13 








159 ± 29 


j 


21 


95±0.05 


-22 


29±0 


05 


185 ± 34 


j 


21 


81±0.06 


-22 


27±0 


05 


172 ± 43 


j 


20 


04±0.06 


-19 


57±0 


17 


115 ± 10 


j 


20 


53±0.07 


-21 


36±() 


05 


80 ± 14 


j 


20 


55±0.06 


-21 


05 ±0 


05 


104 ± 8 


j 


20 


96±0.06 


-21 


27±0 


05 


235 ± 133 


j 


20 


41±0.05 


-21 


16±0 


05 


83 ± 10 




20 


56±0.06 


-20 


23±0 


08 


113 ± 17 




21 


O6±0.05 


-21 


91±0 


05 


204 ± 20 


1 


20 


45±0.05 


-21 


42±0 


05 


118 ± 18 


j 


21 


43±0.05 


-22 


03 ±0 


05 


182 ± 32 


j 


21 


18±0.05 


-21 


51±0 


05 


131 ± 11 


j 


21 


66±0.05 


-22 


98±0 


05 


274 ± 43 


j 


21 


55±0.05 


-22 


13±0 


05 


182 ± 11 


j 


20 


43±0.05 


-20 


76±0 


06 


135 ± 27 




21 


36±0.20 








121 ± 11 


j 


22 


18±0.05 


-22 


83±0 


05 


177 ± 6 




22 


69±0.07 


-21 


99±() 


20 


215 ±3 




21 


60±0.05 


-22 


75±0 


05 


180 ± 11 




22 


12±0.05 


-22 


88±() 


05 


187 ± 19 




21 


33±0.05 


-20 


84±() 


08 


142 ± 11 




21 


53±0.05 


-21 


70±0 


05 


209 ± 10 




20 


51 ±0.09 


-20 


69±0 


27 


74 ± 13 




21 


67±0.06 


-22 


66±0 


05 


174 ± 22 




20 


97±0.10 


-21 


55±0 


19 


128 ± 5 




22 


00±0.05 


-22 


16±0 


05 


172 ± 15 




21 


59±0.06 


-22 


20±0 


12 


59 ± 8 




21 


76±0.06 


-22 


80±0 


05 


258 ± 45 




21 


90±0.06 


-22 


83±0 


05 


194 ± 45 




21 


26±0.07 


-22 


08±0 


14 


146 ± 10 




20 


93±0.05 


-20 


97±() 


28 


157 ± 14 
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TABLE 3 — Continued 



Object 


Sample 


z 


RA 


DEC 


r s 


( 


SI PA 


My 


M K 


y MAX 


Q 








[deg] 


[deg] 


[kpc] 


[deg] 


[deg] 


[mag] 


[mag] 


[km/s] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 




N37427 


field 


0.579 


73.516296 


-2.988545 


5.59 


31.0 


43 


-22.45±0.06 


-23.48±0.05 


215 ±53 


l 


p34i2()c3 


field 


0.595 


6.533985 


16.984360 


4.26 


53.3 


2 


-21.15±0.05 


-21.06±0.32 


109 ±7 


l 


N41286 


field 


0.600 


73.466621 


-2.966479 


2.53 


43.0 


22 


-21.33±0.05 


-20.67±0.11 


107 ± 12 


l 


p34i94c2 


field 


0.614 


6.553948 


16.984819 


11.47 


38.2 


30 


-21.23±0.07 


-21.58±0.33 


251 ± 37 


2 



NOTE. — Cluster galaxies arc arranged in order of increasing R/RyiR- where R VIR = 1.70 Mpc for CI 0024 and R vlR - 2.66 for MS 0451. The cluster centers are (6.6500,17.1433) and (73.5454,-3.0186) J2000, for CI 0024 and MS 0451, respectively. 
Field galaxies are arranged in ascending redshift order. Col. (1): Object name. Col (2): Indicates which subsample that each galaxy belongs to: CI 0024, MS 0451, or field. Col (3): Redshifl of each galaxy. Cols (4) and (5): J2 000 coordinates for each 
galaxy. Cols (6) and (7): Scale lengths r s and inclinations ; measured via GALFIT. (" - 90° indicates edge-on. * Denotes objects where Rij was set manually to improve the fit. Col (8): Indicates the misalignment between the PA of the spectroscopic slit 
and the major axis of the galaxy. Typical formal errors on r s , i, and SI PA are 10%, 0.5°, and 1 ° , respectively; systematic uncertainties in i and SI PA are typically ~ 10° . Cols (9) and (10): Absolute magnitudes in rest-frame V and K bands, corrected 
for inclination-dependent internal extinction. Col (11): Measured V^AX- In cuses where more than one emission line was measured, this is a weighted average. Col (12): Rotation curve quality. Q - 1 curves display turnovers on both ends; Q-2 curves 
display only one turnover, or an uncertain fit. 



